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Abstract 

Recently, studies of the gravitational collapse of a scalar field within spherically sym- 
metric AdS spacetimes was presented in [H [2] which showed an instability of pure AdS 
to black hole formation. In particular, the work showed that arbitrarily small initial 
configurations of scalar field evolved through some number of reflections off the AdS 
boundary until a black hole forms. We consider this same system, extended to include 
a complex scalar field, and reproduce this phenomena. We present tests of our nu- 
merical code that demonstrate convergence and consistency. We study the properties 
of the evolution as the scalar pulse becomes more compact examining the asymptotic 
behavior of the scalar field, an observable in the corresponding boundary CFT. We 
demonstrate that such BH formation occurs even when one places a reflecting bound- 
ary at finite radius indicating that the sharpening is a property of gravity in a bounded 
domain, not of AdS itself. We examine how the initial energy is transferred to higher 
frequencies -which leads to black hole formation- and uncover interesting features of 
this transfer. 
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1 Introduction 



Two recent papers [H E] study the gravitational collapse of a real scalar 
field in spherically symmetric Anti-de Sitter (AdS) spacetimes finding a re- 
markable instability to black hole (BH) formation. That is, given any non- 
vanishing initial scalar field, they provide convincing numerical evidence, 
further supported by analytic arguments, that the resulting spacetime forms 
a BH within a time that scales inversely with the square of the initial am- 
plitude of the scalar perturbation to AdS0. The work discussing this very 
interesting result, however, presents relatively few details and results from 
related work [U |5] are in tension with it. We have therefore reproduced 
aspects of [H [2] and provide new insights in this paper. 

Such a result is, at first sight, surprising from a few different perspectives. 
Since Minkowski and de Sitter spacetimes are stable, it might be expected 
that AdS is stable as well |6]. Within the context of black hole critical 
phenomena [71 [8], the result is similarly surprising. In studies of critical 
behavior, one generally evolves a series of pulses with different amplitudes as 
the initial pulses implode through their centers, searching for the threshold 
of black hole formation. However, this instability suggests that a black hole 
is always formed and that this black hole threshold is only a threshold for 
immediate collapse. 

In contrast, one can similarly view this instability from the perspective of 
the corresponding boundary conformal field theory (CFT). AdS-CFT asserts 
that BH solutions in the bulk correspond to thermal states of the CFT. In 
this light, the fact that this instability would always form a black hole would 
simply imply that states thermalize, which is not surprising at all. Further- 
more, evidence for the instability has also been given in [9l [10] supporting this 
behavior based on ergodic arguments or linearized perturbation analysis as 
well as a theorem guaranteeing the instability under certain conditions [TT] . 
On the other hand, recent works indicate that some scenarios never form a 
black hole, thus thermalization within the gauge/gravity context need not 
be a forgone conclusion [121 US] • 

Furthermore, a pure state can never thermalize during a unitary evolu- 
tion, and so the BH formation conjectures of [H [2], as interpreted within 
the framework of AdS/CFT correspondence, indicate that it might be im- 
possible to set-up a pure initial state in strict limit large- gauge theories. 
Likewise, if the initial configuration carries a conserved global charge, one 
expects (at most) thermalization in a given charge sector. A global charge 
in the boundary CFT dynamics is realized through a local gauge symmetry 
of the dual gravitation bulk dynamics. We postpone the study of thermal- 

^It is interesting to note that Pretorius & Choptuik foresaw this possibihty in their 
study in 2 + 1 AAdS [3]. 
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ization of charged configurations in AdS to future work, and discuss here the 
effect of global symmetries in the bulk on the thermalization. To this end we 
extend the work of [H E] to a complex scalar field collapse in AdS. 

This work is organized as follows; formulation of the problem is given in 
Section O We present a number of details of our numerical work in Section [3] 
and provide tests that our code produces convergent and consistent solutions 
through a number of reflections off the boundary. We consider the results of 
our simulations in Section HI extending to a complex scalar field. In Section [5] 
we interpret the results of our simulations in a context of a dual boundary 
conformal gauge theory dynamics. We conclude in Section El 

2 Formulation 

Consider a. {d+ l)-dimensional gravitational system, dual to a strongly cou- 
pled CFTd on {d — l)-dimensional sphere S''^"^ We focus on S'0((i)-invariant 
states of a conformal field theory with nonzero expectation values for a pair 
of exactly marginal operators Ol^\ i = 1,2. We are interested in a uni- 
tary evolution of such states as determined by the following dual effective 
gravitational action 

where = 0i + z 02 is a massless complex scalar field, dual to a pair of 
exactly marginal operators O^^^ + i O^f^ = Od- Further, 

Md+i = dMd+i X X, dMd+i = Rt^ S"-^ , X = {x e [0, |]} . (2.2) 
Notice that the effective action (12.1 p is invariant under a global phase rotation 

e"*" , (2.3) 
corresponding to a phase rotation of expectation values 

{Od) ^ {Od) e-^°. (2.4) 
Associated with (12. Sp there is a conserved current 

= tg^- (09.0* - 0*9,0) , (2.5) 

and a corresponding conserved charge 
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The other conserved quantity of a CFT during the evolution is the total 
mass M. We can rephrase the work [H [2] as a statement that nonequilibrium 
states {M, Q = 0} of a strongly coupled CFT with a gravitational dual 
evolve into thermal states. Here, in addition to reproducing these claims 
(and reinterpreting the gravitational bulk dynamics from the boundary CFT 
perspective), we are interested in a broader question: what is the role of a 
non-zero charge Q in the process of thermal equilibration? Such a question 
is particularly interesting as it is well-known that a complex scalar field in 
AdS supports perturbatively-stable boson star solutions [T¥j . 

To proceed, we adopt a notation similar to that of |2] choosing the same 
form of the {d + l)-dimensional metric describing an asymptotically AdS 
spacetime 

ds^ = (-Ae-^'de + ^ + sin^ x dnl_^ , (2.7) 

cos X \ -Ti / 

where £ is the scale-size of the AdS spacetime, dVl?^_i is the metric of S"^~\ 
and A{x,t) and 5(x,t) are scalar functions describing the metric. The scale 
^ drops out of the resulting equations and we set it to unity without loss of 
generality. The spatial coordinate x runs over — )■ -7r/2 such that spheres 
through the point x have a radius r = tanx. To describe the real and 
imaginary components of a scalar field, we introduce the two quantities $j = 
{d/dx) {(pi) and = A-^e^ (d/dt) (0^). 

Besides the extension from a real scalar to a complex scalar field, this 
is precisely as described in |2] and we derive the same equations of motion. 
However, working with these variables, achieving a stable numerical scheme 
required a few numerical tweaks at the outer boundary to provide convergent 
evolutions even through reflections (or bounces). A more elegant solution, 
and one which requires no such tricks, involves rescaling the matter functions 
in keeping with their appropriate boundary fall-off. In particular, we use 
rescaled quantities according to 

(2.8) 



cos"' ^ X 



Hi = — — -p:^ = T^, 2.9 

A cos''"^ X cos""^^ X 



cos'^ ^ X cos"^ ^ X 



In terms of these rescaled quantities (we drop the caret from here for- 
ward), we have the following evolution equations resulting from the Klein- 
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Gordon equation 

0i =Ae-^Ui , 

$i = — (cos"^"^ xAe-^Ui) 

sm a; V cosx 
The system also includes two spatial equations to be integrated 



(2.11) 



A,^ = : (1 - ^) - sinxcos^'^-^xA — ^ + n 



5 ^ = - sin a; cos^"*^^ x i — i h 11,^ 



cos^ X 



(2.12) 



together with one constraint equation 

A^t + 2 sin X cos'^'^-^ A^e-^ i^iUi) = 0, (2.13) 
where a sum over 2 = {1, 2} is implied. 

At the origin, these quantities behave independently of d as 

Mt,x) =(pl;\t) + o{x^), 

A{t,x) =l + 0{x^), (2.14) 
6{t,x) =5o{t) + Oix^) . 

At the outer boundary x = 7r/2 we introduce p = 7t/2 — x so that we have 

0.(t,p) =0«(t)p + O(p3), 

^(t,p)=l-M^^ + 0(p2'^), (2.15) 
cos"^ ^ p 

5(t,p) =0 + O{p"'). 

We note here a difference between our work and [2]. We choose conditions 
on 6 such that it becomes zero on the boundary and hence the coordinate 
time with which we examine the boundary describes proper time there as 
well. The gauge freedom in this system allows for such a rescaling without 
affecting the dynamics observed. 

The asymptotic behavior fl2.15p determines the boundary CFT observ- 
ables: the expectation values of the stress-energy tensor T^i, and the opera- 
tors of 

a-1 (2.16) 

16nG,+i{Of) = M 4>f{t), 
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where Ed-, up to an additive constantn for even rf, is proportional to M; g^p 
is a metric on a round S'^~^. Explicitly, for (i = 3, 4 we have 



M, 

3(1 + 4M) 



d = 3 , 
d = A. 



(2.17) 



Additionally note that the conserved U{1) charge is given by 

I-k/2 

Jo 

(2.18) 

Note that since dtQ = 0, the integral in fl2.18p can be evaluated at t = 0. 

The constraint f l2.13p implies that M in fl2.15p is time-independent, en- 
suring energy conservation 

dt£d = 0- (2.19) 
It is convenient to introduce the mass aspect function (t, x) as 



A{t, x) = 1- M{t, x) ^"^_f 

sin X 



Following (12:T21) we find 
M{t,x)= [ dz tan-^-i z cos2('^-^) z A(t, z) 



COS' z 



Comparing (I2.2ip and fl2.15p we see that 

M = M{t,x) 



(2.20) 



. (2.21) 



(2.22) 



Both Sd and the charge Q are determined from the initial data. Motivated 
by P, in this paper we consider the following set of initial conditions: 
■ Q = case 



$^(0,x) =0, 



Ui{0,x) = -e" 

71 



COS 



1-d 



X 6l , 



(2.23) 



where here 61 represents the Kronecker delta. Numerically solving the spatial 
ODEs of ([2J0D and ([MID, we compute A^(0, 7r/2) = M. For c/ = 3 ■ ■ ■ 6 and 
o' = jq the values M(e) are presented in Fig. [H 
■ Q 7^ case 



$i(0,x) = 9, 



— e"^^ cos^-'^x (5. 

TT 



2e 



(2.24) 



nj(0,x) = LOd — e 'r2<T2 cos X 6^ , ujd = d . 

TT 



^This constant determines the Casimir energy of the CFTrf. 
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Figure 1: Mass M of initial configurations f l2.23p as a function of e for 
d = 3, 4, 5, 6 (blue/red/green/orange curves). Note that M oc for Ine < 4. 
While the coordinate extent of the scalar profile is the same in all d, higher 
dimensional AdSd+i more efficiently "localizes" it, resulting in smaller Md 
for a fixed e. 

where Ud is the lowest frequency of the linearized scalar perturbations of 
AdSd+i. Numerically solving the spatial ODEs of (12120]) and dMlD, we 
compute A^(0, 7r/2) = M. For (i = 3 ■ ■ ■ 6 and a = , the values M(e) and 
^ are presented in Fig. [2l 

Note that for both classes of initial conditions, ( 12.23^ and ( ]2.24p . the 
values 0^*''(O), and correspondingly the expectation values of the boundary 
CFT operators O^^ at t = 0, are zero. 



3 Numerics 

3.1 Implementation 

To solve these equations, we employ finite-difference approximations within a 
method-of-lines (MOL) approach that employs a third-order accurate Runge- 
Kutte (RK3) time evolution scheme. Such an approach computes the spatial 
derivatives (with a high-order accurate, finite-difference stencil) and then 
considers the evolution equations as ordinary-differential equations (ODEs) 
in the time coordinate. 

As is standard in studies of critical behavior, we use adaptive mesh re- 
finement (AMR) to add resolution where needed and remove it when not 
needed. In particular, we use Choptuik's ad infrastructure [3 US] adapted 
for our uses here via the changes: 

• modified the loop for convergence within the iterative Crank-Nicholson 
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Figure 2: Mass M of initial configurations ( I2.24p as a function of e for 
d = 3 ■ ■ ■ 6 (left panel, blue- ■ ■ orange) and corresponding ratios ^ (right 
panel). Here, M oc scaling of mass extends only up to In e ~ 1. Note that 
^ approaches a constant as M — > 0, while In ^ oc M for large M, effectively 
making the charge negligible. The rescaling of the horizontal axis of the right 
panel by 5*^"^ is made for readability of the graph. 

scheme to a three-step RK3 update 

• replaced second-order accurate spatial derivatives with high-order ac- 
curacy spatial derivatives 

• modified the AMR-boundary treatment from interpolation for a single 
point in time to interpolation in time of width two points in the parent, 
with spatial interpolation on the fine level between these two points. 
This approach follows in spirit the tapered approach of [16], although 
it strictly lacks the higher-order convergence there. 

The numerical grid has two boundaries, the origin and the AdS boundary 
at p = 7r/2 — X = 0. To enforce the boundary conditions (12.141) at the origin, 
we set $j(t, 0) = and A{t, 0) = 1 and obtain values for (f)i{t, 0) and Ili{t, 0) 
from quadratic fits. Because 6{t, x) is computed by numerical integration 
from the outer boundary, no boundary condition is needed at the origin. At 
the p = boundary, we enforce the conditions (I2.15P by setting 0j(t,7r/2), 
Ili{t,7r/2), and $j(t,7r/2) all to zero there. We set the metric functions 
A{t, 7r/2) = 1 and 6{t, 7r/2) = 0, the latter chosen to ensure that coordinate 
time t measures proper time on the boundary. 

3.2 Tests 

We carry-out a number of tests of our implementation. We find that the 
solutions obtained generally converge at least as fast as third-order in the 
grid spacing. We evaluate the residual of the momentum constraint and 
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confirm that it approaches zero as the resolution is increased. We compute 
both the total mass A^(t,7r/2) using (12.211) and charge using fl2.18p . and 
ensure that their degree of conservation in time improves with resolution. 
We also completed two standard black hole threshold searches and find the 
expected Choptuik behavior [7j. Fig. [3] presents some results from one such 
test. For the AMR evolutions mentioned here, we ensure that increasing 
both the coarse level resolution and decreasing the threshold for refinement 
produce essentially unchanged collapse times. 

4 Results 

We find the same behavior as found in [Tl[2]. In particular, given any initial, 
non-vanishing configuration of scalar field, an evolution with sufficient reso- 
lution and time produces a black hole. We observe that as the pulse travels 
back and forth from the boundary, one sees increasing compactness in the 
scalar profile. Examining the metric function A{x, t), one sees that on return 
of the pulse to the centre, the metric function approaches to an increasing 
degree, and it is precisely the approach of A{x, t) to zero that indicates the 
formation of a black hole. In other words, the increasingly compact scalar 
profile produces an increasingly deeper gravitational potential well. The lo- 
cation of this minimum also decreases with each bounce. We plot these two 
quantities in Fig. IHfor one particular run of the initial data of f l2.23p . 

4.1 The Nature of the Instability 

One can ask whether this instability results from the particular nature of 
AdS spacetime or instead just the fact that one is evolving in an essentially 
bounded domain. Beyond the perturbative analysis, one can imagine sending 
a scalar pulse down an elevator shaft equipped with a mirror. It will blueshift 
as it gets deeper, but upon return to its original height, it would obtain its 
original form only to first order. At higher order, one considers not just the 
gravity of the Earth, but also the self-gravity of the pulse which will act to 
compress the pulse to its center. 

One could confirm this picture by numerically evolving a scalar pulse in 
an asymptotically flat space but with a reflecting boundary at some radius, 
as done in [17]. Instead, we do something similar but with our same code. 
That is, we enforce reflecting boundary conditions on the scalar field at some 
a^finitc < 7r/2 

0i(a;finitc, ^) = , nj(Xfinitc, t) = 0. (4.1) 

The value is arbitrary, but by truncating the domain, the scalar pulse does 
not "see" that the spacetime in which it lives is AdS. We choose ajgnite = 
7r/4 and find that pulses behave qualitatively in the same way, i.e. as they 
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Figure 3: Convergence test for two bounces with d = 3. The initial data is of 
the form (12.231) with e = 20.01 and a = 1/16. (top) The order of convergence 
obtained from comparisons of the metric function A{x, t) at successively dou- 
bled resolutions. The convergence order is computed from three different 
resolutions, a base resolution and runs with half and one-quarter the base 
grid solution. Thus, run "32" has a grid spacing 2~^ that of the run labeled 
"1." (middle) Order of convergence obtained for <I>i(x, t). Both these results 
indicate convergence is better than third order convergent, (bottom) The 
logarithm of the L2-norm of the momentum constraint residual for just the 
three best resolutions. That it decreases with increasing resolution suggests 
the results are converging to a proper solution of Einstein's equations. Note 
that this residual is computed only at first order of accuracy because it in- 
volves a time derivative and therefore one should not estimate the order of 
convergence from it. 
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bounce, they shift to shorter wavelengths and eventually form black holes. 
We subsequently varied the value of Xfinitc and confirmed that this behavior is 
robust. Interestingly when e in initial data of the form (12.23^ was decreased 
below 18.39, collapse was not observed even after hundreds of bounces. 

4.2 Boundary Information 

Given the AdS-CFT correspondence, it is interesting to consider what the 
boundary CFT sees from the bulk. One usually uses a "dictionary" to ex- 
tract various quantities. However, here the boundary stress-energy tensor is 
time-independent and just equal to the total mass of the spacetime, M. Of 
particular relevance is the leading asymptotic value of the scalar field 0j as its 
behavior corresponds to a quantum operator on the boundary. We therefore 
determine the asymptotic behavior of each scalar field 0j with a polynomial, 
least-squares fit to 

Up.t) = ct>f{t)p + <P%,{t)p\ (4.2) 

As a check of our extracted value, we employ from ten to thirty points and 
confirm the results of the relevant quantity, cf)"^ (t) are hardly affected. 

This asymptotic information is displayed for a particular case in Fig. |5l 
The left side shows the behavior as a function of time with each successive 
bounce along with the Fourier power spectrum. This spectrum is computed 
usmg the FFTW library [18] acting on (j)^^\t) once it has been interpolated 
onto a uniform grid over t (recall the use of AMR produces nonuniform in 
time information for the asymptotic behavior). The power is computed as 
simply the square of the (complex) Fourier amplitude. 



5 Boundary CFT perspective of the gravitational col- 
lapse 

The formation of an apparent horizon (AH) in AdS gravitational collapse 
is holographically dual to the thermalization of an initial CFT state cor- 
responding to the initial condition of the gravitational evolution. To fully 
characterize the gravitational dynamics from the boundary perspective, one 
would have to compute an infinite set of local correlation function^. Here, 
we focus on the simplest local observables, namely (O^*^) oc Recall that 
the one-point function of the boundary stress-energy tensor is constant dur- 
ing the evolution (see ( 12.16^ ). and thus does not carry any information about 
AH formation. The same applies to the conserved U (1) charge and so its role 
is to identify (select) dynamical thermalization trajectories. 

■^Also, presumably, the nonlocal observables such as Wilson loops. 
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Figure 4: Demonstration of the increasingly compact gravitational poten- 
tial with successive bounces of initial data of form f l2.23p with e = 20 and 
(T = 1/16. Six bounces are shown before a black hole forms, (top) The 
minimum of A{x, t) as a function of time. At the final time, the minimum 
approaches zero, signaling black hole formation, (middle) The a;-coordinate 
where A{x, t) achieves its minimum. In these coordinates, the speed of light 
is unity, and the line segments roughly indicate the motion of the scalar pulse 
back and forth across the grid. That compaction increases with each bounce 
is quite hard to observe on the linear scale, and so the logarithm of this data 
is shown at bottom. The evolution terminates at the last time shown as the 
solution approaches black hole formation. 
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Figure 5: The behavior of the scalar field at the boundary with successive 
bounces for the same evolution as in Fig. HI At left, the asymptotic value 
(f)^^\t) from Eq. (14.21) is shown for segments of time of length vr. The pulse 
sharpens with each "implosion" through the origin. On the right is shown 
the Fourier transform of the signal on the left. Shown is the natural logarithm 
of the amplitude of the Fourier component squared as a function of (the log 
of) angular frequency u. The vertical scale is arbitrary but fixed for all 
bounces. The transform shows that energy is shifted to higher frequencies. 
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5.1 Weakly nonlinear perturbations with global charge in AdS 

Following P (see also [6]) we consider the solution of ( ]2.1ip - (l2.13p . pertur- 
bative in the bulk scalar amplitudes e 



oo 



^ e^^+i , A = 1 - ^ e^^- ^2,- , b = Y, e'^ 82, , (5.1) 

i=o j=i 3=1 

where 0i,2j+i , 5(^2^ are functions of (t,x). It is convenient to decompose 
these functions in terms of a complete basis. A natural basis is provided by 
the AdSd+i massless scalar eigenvalues and eigenfunctions (which we refer to 
from now on as oscillons) 

ujj = d + 2j , ej{x) = dj cos'^ x 2F1 (^-j, d + j, ^, sin^ , j = 0, 1, ■ ■ ■ , 

(5.2) 

where dj are normalization constants such that 

7r/2 

dx ei{x)ej{x) tan*^""^ x = 5ij . (5.3) 





A remarkable observation of [1] was that initial conditions which represent 
at a linearized level (at order 0{e) ) a superposition of several oscillons with 
different index j appear to be unstable at time scales tinstabiuty ~ C'(e~^); on 
the other hand, nonlinear effects of a single oscillon do not lead to destabi- 
lization. Specifically, the instabilities occur whenever oscillons with indice^ 
{ji)j2,j3} are present at order 0{e), while the oscillon with index jr, such 
that 

is not excited at this order. 

In what follows we argue that the basic conclusion of [1] remains valid 
even for configurations with Q ^ 0. We point out though an interesting 
twist: although any finite number of distinct-index oscillons excited at order 
0{e) leads to nonlinear instabilities at order 0{e^), 0{e) initial configurations 
with an infinite number of distinct-index oscillons remain (formally) stable 
at order O(e^). 

While it is straightforward to carry out the analysis in full generality, to 
be maximally explicit we discuss the case of c? = 3 and a superposition of two 
lowest-frequency oscillons at a linearized level, namely Cq and Ci. We discuss 
different cases: 

• a single neutral oscillon 



*The indices could be repeated. 
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• a superposition of two neutral oscillons 

• a single charged oscillon 

• a superposition of two charged oscillons 

• a superposition of a charged and a neutral oscillon. 

Lastly, we comment on nonlinear instabilities at 0{e^) arising from different 
channels (15.41) . i.e., different combinations of {ji, j2) Js} resulting in the same 

To distinguish neutral and charged initial conditions it is convenient to 
express the charge Q given by fl2.18p as 



g = 47rM dx tan2a;A"^(t,a;) e^(*'^) ((9j0i(t, a;)02(t, - 9t02(t, a;)0i(t, a;)) . 



Mt,x)=Mt,x), 0i(O,x) = eeo(x) + 0(e=^), 9t0i(O, x) = . (5.6) 



Clearly, in this case Q = 0. Requiring normalizability and regularity of 0i, 
to order 0{e^) (where first nonlinear effects appear), we find 




(5.5) 



5.1.1 A single neutral oscillon 



Here we take 




135 
An 




(5.7) 



+ F3,9(x) cos(9t) +C(e^), 



with 



3,3 




12 COS® X — 88 cos^ X + 108 cos"^ a; — 63 cos^ x + 637r^ 




(5.8) 



Fa^g = , cos^ x(9 cos^ x — 4) . 



(5.9) 



Notice that in (15. 7p we absorbed a term linearly growing in time 



oc e^t sin(a;ot) 
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into O(e^) shift of the leading-order oscillon frequency Uq: 

uo^wo-^e"^. (5.10) 

Obviously, we could do so because an oscillon with such a frequency has al- 
ready been present in the initial condition (15. 6p . For this initial configuration 
the instability condition (15. 4p is satisfied only for ji = j2 = js = jr = 0- 



5.1.2 A superposition of two neutral oscillons 

Consider now a slightly more general (neutral) initial condition 

02(t,x) = 0i(t,x), 0i(O,x) = e(eo(x)+ei(x))+0(e3), 9t0i(O,x) = O 

(5.11) 

Once again, in this case (5 = 0. Requiring normalizability and regularity of 
01, to order C(e^), we find 



, , o 335 2\ A , . ( ( ^ 1519 2 , 
eo(a;) cos ( ( 3 - — e 1 M + ei(a:) cos 115 —e \ t 



+ 



8 

^F3,2fc-i(x) cos((2A: - l)t) + ^e2(x) t sin(7t) 



fc=i 



0(6^) , 



(5.12) 

where ^3^2^+1(2^) are some analytically determined function^. Here, we have 
three different terms at order 0{e^), which grow linearly with time 



oc e^t X < cos(a;ot) , cos(a;it) , sin(ci;2t) 



(5.13) 



The presence of j = {0, 1} oscillons in order 0{e) initial conditions allows 
us to absorb the first two terms into the shifts of the leading-order oscillon 
frequencies 



Uo Wo 



335 



1519 



, ui^wi- -^e^. (5.14) 
Sir OTT 

We cannot do the same with the remaining term in (I5.13P — for this to 

happen 0i(O,a;) must contain a term oc e 62(2;). Of course, the presence of 

62(3;) at order 0{e) in the initial conditions, while eliminating e^t x sin(ci;2t) 

term, would generate new resonances at j > 2. Thus, the absence of linearly 

growing with time terms at order 0{e^) necessitates the excitation of all 

oscillons at order 0{e) in the initial condition. In what follows we show 

that 0{e^) instabilities can be removed by shifting the leading-order oscillon 

frequency (provided the relevant excitation is present at order 0{e)), even if 

they arise from different channels (15. 4p . 



^We omit their explicit expressions to keep formulae readable. 
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5.1.3 A single charged oscillon 
For charged initial data we take 

0i(O, x) = eeo(x) + 0{e^) , 9i0i(O, x) = , 
02(0, = , 9t02(O, x) = eujo eo{x) + ^(e^) 

where (see (15. 5p ) 

Q = -127r2e2 + C(e^) . 



(5.15) 



(5.16) 



Initial data f l5.15p is a direct charged generalization of a single-oscillon initial 
data f l5.6p . In this case, to order 0{e^), we find 



h = e 
^2 = e 



eo{x) cos ( ( 3 - ) t 
eo(x)sin ( ( 3 - ) t 



Fl,{x) cos(3t) 
F|,(x) sin(3t) 



C(e' 



0(6^) . 



(5.17) 



All the linearly with time growing terms are absorbed into the higher-order 
shift of Uo- 



{5.11 



5.1.4 A superposition of two charged oscillons 
For charged initial data we take 

0i(O,x) = e(eo(x) + ei(x)) + 0(e=^) , 9t0i(O,x) = 0, 
02(0, x) = , 9j02(O, x) = e{ujQ eo(x) + Ui ei(x)) + C(e=^) , 

where (see (15. 5p ) 

Q = -^2tx\^ + 0{e^) . (5.19) 

Initial data (I5.18P is a direct charged generalization of the two-oscillon initial 
data (15. lip . In this case, to order 0{e^), we find 



>i = e 



159 



2201 



eo(x)cos^^3 e J tj + ei(x) cos ^^5 ^ 

^ 50\/6 
5^^2^-1(3;) cos((2A:- l)t) + ^e2(x) t sin(7t) 



0(e^ 



!>2 = e 



eo(x) sin I 13 



159 



2201 



r 4 



E^3Vi(^) sin((2fc-l)t) 



fc=i 



50^6 
37r 



62 (x) t cos(7t) 



0(e^ 



(5.20) 
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Parallel to the discussion in Section I5.1.2[ the absence of leading order oc 
ee2{x) modes in the initial condition flS.lSp leads to growing resonance terms 



oc e^t X |cos(a;2t) , sm{u2t) 

5.1.5 A superposition of a charged oscillon with a neutral one 
For charged initial data we take 

0i(O, x) = e(eo(x) + ei(x)) + ^(e^) , d^Q, x) = , 
02(0, x) = , 9t02(O, x) = euo eo(x) + ^(e^) , 

where (see (15. 5p ) 

Q = -Unh^ + C(e^) . 
In this case, to order 0{e^) we find 



(5.21) 



(5.22) 



(5.23) 



h = e 



787 



1843 



e'\t 



eo(x) cos ( ( 3 - iT^e^J t j + ei(a;) cos (^(^5 

8 or fF 

Yl ^3Vi(^) cos((2A; - l)t) + ——e^ix) t sin(7t) 



?2 = e 



eo{x) sm 



153 



-e'\t 



E ^3Vi(^) sin((2A; - l)t) 



75 5^/6 

—ei[xj t cos(5t) + 62 (x) t cos(7t) 

87r 127r 



0(e^ 



(5.24) 



Parallel to the discussion in Section I5.1.2[ the absence of leading order oc 
ee2{x) modes in the initial condition f l5.22p for both (pi and 02 leads to growing 
resonance terms 

oc e^t X < cos(w2^) , sin(cj2^) \ ■ (5.25) 



Likewise, the absence of a leading order oc eei(x) mode in the initial condi- 
tion f l5.22p for 02 leads to a growing resonance term 



oc eH cos(a;it) . 



(5.26) 



5.1.6 Multi-channel instabilities at 0{e^) and their removal 

Previously, we considered the linear (order 0{e)) superposition of two neu- 
tral/charged oscillons and identified instabilities at order 0{e^) arising from 
the resonance term with index jV (see f l5.4p ). Furthermore, we showed that if 
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an oscillon ej^{x) is present in the initial data at order 0{e), these instabihties 
can be removed with an appropriate 0{e^) shift of the oscillon frequency 

u,^ -> u,^ + . (5.27) 

Thus, it appears that to avoid 0{e^) instabilities, initial data at order 0{e) 
must contain all oscillons. 

In all examples discussed, only a single instability channel was present. 
If all oscillons are excited at order 0{e), the same resonance frequency Wj^ 
(15. 4p would arise from different channels. We argue here that all such multi- 
channel instabilities can still be removed with a single resonance frequency 
shift (I5.27p . To be maximally explicit, consider a neutral, order 0{e), initial 
data containing oscillons j = {0, 1, 2, 3} 

3 

02(t,x) =0i(t,x), 0i(O,x) = e ^Aje,(x) + C(e3), dtMO,x) = 0, 

j=0 

(5.28) 

where Aj are generically different amplitudes. Consider a resonance index 
jr = 2. From (15. 4p there are 7 distinct instability channels: 



(1): 


Wo 

Jr 


= U2- 


h Wo 




(2): 


OOi 


= UJ2+UJ1 


- UJi 


(3): 




= 002- 


h 002 




(4): 


^jr 


= OJ2+OJZ 


- 


(5): 


^Jr 


= Ui - 


h OOi 




(6): 


^jr 


= + Wo 


- UOi 


(7): 


^jr 


= UJ3- 


h OOi 


- 0O2 . 











(5.29) 



Each of the channels in (I5.29P would contribute a term that grows linearly 
with time at order 0{e^) to (pi 

c iro^r,u„uro ■?,-/, N / N 1 /60613 .o 35^/6 . .n 189\/5 . . . 

= eHsm{uj2t) e2{x) -( + ^^^3^5 + -^A^AgA^ 

50607 . . 1090 ,2 . 1797 ,2 . 159^30 , , , 

H AiAr H AiAr H AiAj H A^AjAc 

40 ^ 3 ^ ^ 20 2 ^ ' - 



(5.30) 



A resonance frequency shift (15.27P would additionally contribute 



^^shift _ eHsm{oo2t) e2{x) ^(^Ar^^^j • (5-31) 
Clearly, provided Aj 7^ 0, all the instabilities (I5.29P can be removed, i.e. 

^^rcsonancc ^ ^^shift ^ q ^ 

for suitably adjusted ^2^^ 
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Straightforward, but quite tedious, considerations show that for a generic 
neutral initial condition 



j=0 

(5.32) 

the presence of resonance contributions from channels 

^jr = + - ' (5-33) 
can be eliminated by the following resonance frequency shift 

^> ^ ^> + E ^^^i^ X 0(1) , (5.34) 

where summation occurs over all possible instability channels (triplets {ji, j2, Js} 
satisfying (15.331) ). The amplitudes Aj are expected to decay fast enough as 
j — oo for the sum in (15.341) to be convergenl0 (and also not to have an 
apparent horizon already present at the initial hyper surf ace). 

To reiterate, the reason why renormalization of the basic resonance fre- 
quencies would remove all instabilities up to order Oie'^ is because to this 
order the necessary frequency shifts due to distinct instability channels (I5.33P 
add-up linearly, as explicitly demonstrated by example (I5.28p . From the lat- 
ter example it is also clear that the 'strength' of each of such shifts — iJ^2^ 
in equation (15.311) — is inversely proportional to the amplitude of a resonance 
frequency at order C(e), i.e., oc -j-. Thus, provided all resonances are excited 

at 0{e) (so that this 'inverse amplitude strength' is finite), the instabilities 
can always be eliminated as (schematically) indicated by (I5.34p . 

Note that while we discussed neutral initial conditions, similar consider- 
ations also apply for initial conditions carrying global charge. Whether or 
not weakly-nonlinear instabilities can be removed (or removed under some 
conditions) beyond order 0{e'^ is an interesting open question. We empha- 
size that our numerical simulations, as well as those reported in [H 12] are 
outside the weakly nonlinear regime. 

5.2 Global charge and the thermalization time 

In global AdSd, the time it takes for a null geodesic to cross the space-time 
is 

^crossing ^ • (5.35) 



^In [13] we present an explicit example where this is the case. 
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Figure 6: AH formation time ^ah a function of mass M for d = 3 (left 
panel) and d = 4 (right panel) on a log- log plot. The blue dots correspond 
to thermalization of states with Q = 0; the red dots correspond to thermal- 
ization of states with Q ^ 0. In both cases, for small M (or e), we find 
that AH is formed earlier than oc (or oc e"^) as suggested by weakly 

nonlinear analysis. A global charge Q of the initial configurations f l2.24p 
does not extend the thermalization time, compared to the same-mass Q = 
configurations (12.231) . 

In Fig. [6] we present the ratio of the AH formation time ^ah to tcmssing as 
a function of mass M. The blue dots represent collapse from initial con- 
dition f l2.23p . and the red dots correspond to collapse from initial condi- 
tion f l2.24p . We find that an AH is forming earlier for Q 7^ initial conditions 
with the same mass as the corresponding Q = initial conditions. This dif- 
ference in collapse times could potentially be attributable to either: (i) that 
one has a global charge while the other is neutral, or (ii) that the spectral 
content of the two initial conditions are different. To understand which of 
these is responsible, we decompose the mass aspect function fl2.20p into the 
oscillon basis fl5.2p 

00 

M{t = 0,x) = Y^ rrij ej{x) . (5.36) 

j=0 

The spectral content for the two initial conditions (both with the same mass) 
is presented in Fig. [71 The spectrum of the charged configuration does not 
show a noticeably broader bandwidth or enhanced amplitudes at high fre- 
quencies when compared to the spectrum of the neutral configuration. Thus, 
it appears that the acceleration of the collapse is a dynamical effect related 
to the class of initial conditions fl2.24p . 

In a similar fashion, we find significant differences in collapse times when 
changing the spacetime dimensions d. However, in this case the spectral 
content of the initial mass aspect function does account for the difference in 
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Figure 7: Decomposition of the mass aspect function A^(t = 0,x) into 
the oscillon basis for two contrasting initial configurations: for e = 10 
neutral fl2.23p initial condition (blue dots), and the charged initial condi- 
tion (12 .24^ with the same mass (red dots). Left panel: d = 3, right panel: 
d = 4. The different spectra in the oscillon basis do not appear to account 
for the difference in the observed AH formation times. 

AH formation time. The left panel in Fig. [S] compares AH formation times 
for d = 3 (blue dots) and d = 4 (red dots). The spectral content of the 
corresponding mass aspect function is shown in the right panel. Because 
higher dimensional oscillons are more localized in the AdSd+i center with 
increasing d, there is larger overlap of the initial profile fl2.23p with higher 
modes for larger d, resulting in broadening of the spectral bandwidth. 

5.3 Bulk vs boundary turbulence 

A motivating question behind this study wonders what the CFT sees of this 
instability on the boundary. Because the only observables on the boundary 
we consider are the (unchanging) total mass M and the asymptotic behavior 
of the scalar fields 4>^^\t) (see (12.151) ) (or equivalently (C^*^), see (12.161) ). we 
examine the spectral content of the latter. 

Fig. E] displays snapshots of the scalar field for successive bounces dur- 
ing an evolution. The initial pulse becomes increasingly more compact, ap- 
proaching something like a Kronecker delta function. As it does so, the FFT 
of each pulse similarly approaches the FFT of a delta function, in particular, 
a flat, plateau-like spectrum. In general, if we define the effective bandwidth 
of the power spectrum as the frequency at which the plateau ceases, then 
we observe this bandwidth to increase rapidly during an evolution. The only 
exception appears, for some runs, during the last few bounces, during which 
the bandwidth decreases modestly from its maximum. However, we specu- 
late that this behavior can be attributed to the nonlinearities of gravity near 
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Figure 8: The AH formation time noticeably decreases with increasing 
boundary space-time dimension d (left panel). Blue dots represent d = 3 
and red dots d = 4. This decrease can be explained by the broadening of 
the oscillon decomposition spectrum (see fl5.36p ) of the mass aspect function 
corresponding to initial conditions for different d (right panel). 

black hole formation. 

Before discussing the behavior of the bandwidth, we should comment on 
the falloff behavior of the power spectrum. It is customary in studies of tur- 
bulence in fluid dynamics to observe the rate of falloff (the slope on a log-log 
plot). In particular, a famous result of Kolmogorov expects a slope in the 
energy spectra of —5/3 within the inertial regime of turbulent flows (inde- 
pendent of spatial dimension d). Here, one might consider the scalar field, 
in analogy with a fluid, as a "gas" of oscillons in which structure transitions 
to increasingly higher frequencies. An analysis of the energy spectra over 
the oscillon basis in the bulk in [T7] finds a fall-off at the rate of —1.22. As 
well, an analysis of the Fourier spectrum over the time series of the Ricci 
scalar at the origin in [19] reports a slope of —1.61, seemingly close to that 
of Kolmogorov. However, the Ricci scalar is an energy density and thus the 
expected scaling is not —5/3. Indeed, we can retrace Kolmogorov's argument 
for such quantity by the following dimensional considerations. The Fourier 
transform of H^ (which scales as £^ oc T~^) scales as £^ oc LT~^ for a typical 
wavelength L. If turbulent behavior transfers energy to different scales, for 
H^ it does so with an efficiency rj which has dimensions of (transfer per 
unit time). Thus, one would expect these quantities to be related by 

(5.37) 

where p and a are (so far) arbitrary powers. However, for fl5.37p to be 
dimensionally consistent, then we must have p = 2/3 and and a = — 1. 

In evaluating what our data says about the falloff, we have to know which 
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parts of the spectrum represent true physics as opposed to numerical noise. 
Straightforward unigrid convergence tests suggests that we can roughly trust 
the spectra up to a bit better than u 500 (i.e. the Nyquist frequency 
un = it/ {2 Ax) with Ax the grid spacing). However, in general, we run with 
AMR and can resolve locally enormously better than any of our unigrid evo- 
lutions. One issue is that while AMR is very efficient locally, its performance 
on global measures is difficult to evaluate. Instead, we carryout some runs 
with different AMR parameters and evaluate where the spectrum changes 
significantly. 

Additionally, AMR introduces its own noise into the asymptotic behavior 
of the scalar field, especially for higher dimensions {d > 6). This noise 
appears as very small kinks in the time series data which contribute noise 
that behaves as oc o;"^ (such that the power fall-off from this noise source is 
-4). 

We also note two other sources of noise in the spectrum. In order to 
carryout a Fourier decomposition, we first interpolate to a grid uniform in 
t. Comparing FFTs with either linear or cubic interpolation, we find this 
noise is generally small, only significant at quite high frequencies and not 
strongly dependent on the interpolation order. The final source of noise 
arises from potential mismatches in the time-series data between (j)~^{mi) 
and ip'd'iV'' + 1] Because the FFT expects a periodic signal, this mismatch 
appears as a (quite small) discontinuity producing noise that falls off as oc uj~^ 
(with a power that falls off with exponent —2). 

As a practical matter, we find that for different runs, different noise 
sources dominate the high-frequency behavior of the power spectrum. Indeed 
the only robust features of the various power spectra we have computed at 
the AdS boundary have been the delta-function plateau and high-frequency 
noise the power of which falls off either with slope —2 or —4. In particular, to 
make a convincing argument that some particular fall-off is characteristic of 
this weakly nonlinear instability, one needs to demonstrate not only a feature 
of the spectrum independent of any noise, but also a feature not present in 
the initial spectrum. That is, any characteristic feature must develop from 
the evolution. 

At the boundary, the only such feature we find is that of the plateau 
and we study its defining feature, its bandwidth. In Fig. [HI we display the 
results of measuring the bandwidth for one set of evolutions. We show the 
power spectrum of the first, middle, and last bounces. From each bounce, we 
calculate a maximum frequency as the location where the spectral power falls 
some threshold below the maximum (typically 20% below the maximum). 
This frequency then serves as an estimate of the bandwidth of the pulse and 
is shown on the right side of the figure. As the pulse sharpens in time, its 
spectrum reaches into higher frequencies and the rate at which it does so is 
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roughly exponential. 

Looking at the right side of Fig. |9l one notices that the growthrate of 
the bandwidth, after an initial quiescent period, appears to approach a line, 
representing exponential growth of the bandwidth. We are thus led to ex- 
amine the bandwidth of runs with a relatively large number of bounces, and 
indeed we find exponential growth of the bandwidth for these runs, shown in 
Fig. [101 

What is particularly interesting about Fig. [10] is that these runs demon- 
strate what appears to be a common growth rate with slope 0.04. This 
commonality in slope appears despite the fact that these families span dif- 
ferent dimensions d and different forms of initial data. What they have in 
common is that they all form black holes after between 26 and 39 bounces. 
Note that the growth rate for these evolutions differs from that shown for 
d = 3 in Fig. M (with a slope of roughly 0.07) which collapses after just 21 
bounces. 

Another interesting aspect of these bandwidth-vs-bounce plots is that 
nearly all evolutions studied terminate at roughly logj^o i^^) ~ 2.7. Presum- 
ably this reflects that black hole formation occurs for such high-frequencies. 
We speculate based on these two features, that a bandwidth plot such as 
discussed here can be considered as something of a phase-space picture for 
this system. At the top of the plot, is black hole formation. Low on the 
left side, evolutions not yet consisting of a black hole evolve through some 
quiescent period. Depending on how weak the initial data, they enter an ex- 
ponentially growing bandwidth regime until they are brought to black hole 
formation. Perhaps, the growth rate in this linear regime is set by how low 
they begin, or, equivalently, how many bounces they require to achieve the 
black hole scale. We stress that this is speculation and will require more 
runs to test whether indeed common growth rates occur for disparate evo- 
lutions that share the same number of bounces. Similarly, this phase space 
picture may require some other basis, such as the oscillon basis, over which 
to compute the bandwidth. 

Last, we also monitor the behavior of 11^ at the origin which is related to 
the energy density (and scales in the same way). Numerically, the behavior 
of the scalar fields at the origin produces less noise and hence cleaner power 
spectra. Figure [TT] presents results for d = 3,4,5, and 6 for the last bounce 
before black hole formation. At frequencies below u ^ 800 a plateau behavior 
results from the pulse profile being quite narrow and approaching a delta-like 
behavior. Within the range u G (800, 3000) the slopes measured fall in the 
interval ~ (—1.8, —2.2), consistent with our previous Kolmogorov-type argu- 
ment for the power spectrum of a quantity scaling as energy density (recall 
that the power spectrum is the square of the complex Fourier amplitude and 
hence should fall-off as —2a). 

The relatively wide range of this interval is perhaps not surprising given 
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that we lack sharp knowledge of where the inertial regime lies in the spec- 
trunl^. On one hand, because energy is only "injected" initially, as it shifts to 
higher frequencies it pushes the initial lower bound for the inertial frequency 
higher. On the other hand, at high frequencies one has a "dissipative" scale 
determined by the black hole size. However, because energy is also "spent" in 
curving the spacetime, the location of this dissipative bound is also difficult 
to determine. Nevertheless, inspection of slope breaks in the low and high 
frequency regimes leave an intermediate domain in which the power spectra 
appears to fall-off with a slope in the range mentioned. Pushing the analogy 
with turbulent fluid dynamics even further, one could speculate that devia- 
tions from a slope of —2 are not surprising, but instead that they would arise 
much as intermittency arises in turbulent fluids. With fluids, intermittency 
introduces deviations from the —5/3 slope due to the presence and dynamics 
of vortices. Here perhaps the sharpening pulse as it periodically propagates 
through the origin could leave a similar imprint on the spectrum. 

6 Conclusion 

We have reproduced the work of [H [2] and provided details of our numerics 
and tests. We find that this instability towards higher frequencies occurs even 
in a finite-sized domain with refiecting boundary conditions that would not 
be expected to know that the spacetime is asymptotically AdS. This result, 
along with the similar case in asymptotically Minkowski in the work of [T7] . 
suggests that the instability is a property not of AdS itself, but instead of 
gravity in a bounded domain. 

We generalized the work of [H |2] to complex scalar collapse in AdS to 
address the question of whether AH formation in global AdS from generic 
initial conditions from the bulk perspective (or thermalization trajectories 
from the boundary CFT perspective ) is sensitive to the presence of a con- 
served global charge. For the class of initial conditions we study (f l2.23p and 
f l2.24p ) the answer appears to be "no." Rather, the collapse is affected by the 
spectral decomposition of the initial data in the oscillon basis (see (15. 2p ). We 
would like to stress that global charge can (and does [13]) strongly influence 
AH formation trajectories in the vicinity of stable stationary conflgurations 
in AdS. 

We repeated the weakly nonlinear analysis of the gravitational collapse 
in the oscillon basis, originally presented in [1] . We argued that to the lead- 
ing order in nonlinearities, all the instabilities can be removed, provided all 
oscillons are excited at the linearized level. This conclusion is robust and 

^To more easily identify this regime, one could introduce boundary deformations to 
continuously pump energy from the boundary into the system as done in e.g. |20j 
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Figure 9: Demonstration of the bandwidth growth with each bounce. Shown 
are the resuhs from the evolution of initial data of form f l2.23p with e = 8, 
a = 0.1 for dimensions ci = 3, 4, 5 and 6. On the left is shown the spectrum 
of the first (blue), middle (green), and last (magenta) bounces. The pulses 
increasingly demonstrate a plateau-like region indicative of becoming more 
like a delta function. At frequencies above the plateau, unphysical noise is 
apparent. On the right is shown the bandwidth associated with all pulses. 
The bandwidth is obtained as the location at which the spectrum drops from 
its maximum by twenty percent. As noted in the text (see Fig. [8]), the same 
initial profile for the scalar field contains higher oscillons in its decomposition 
with higher dimensions, and therefore one observes the bandwidth of the 
pulses grow more quickly than the d = 3 case. 
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Figure 10: Bandwidth as a function of bounce for a few different evolutions. 
Choosing evolutions with roughly 30 bounces with different dimension d and 
different families, we observe similar, exponential growth of the bandwidth 
for some intermediate regime (roughly from bounce 10 to bounce 30). Shown 
are the results from: (blue crosses) initial configuration fl2.23p with e = 10 in 
d = 3, (green open circles) initial configuration f l2.24p with e = 10in(i = 3, 
(cyan open circles) initial configuration f l2.23p with e = 7 in c? = 6, (magenta 
triangles) initial configuration (I2.23P with e = 7 in = 5, and (red squares) a 
Gaussian initial profile m d = ?>. four (black dashed) lines all with the slope 
0.04 but with various y-intercepts chosen to align with the bandwidth data 
are shown. 
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Figure 11: Power spectra of 11^ at the origin versus frequency u for d = 
3,4,5,6 as measured in the last bounce before black hole formation. At 
frequencies below ^ 800, the apparent plateau behavior is indicative of the 
sharpened profile of the pulse. At higher frequencies up to w ~ 3000, the 
spectral power decreases with a slope in the range G (—1.8, —2.2). 
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independent of the global charge. We further argued that, while the picture 
of gravitational collapse as an ideal gas of oscillons is valid in the linearized 
regime, this interpretation is in conflict with the late-time evolution. 

Finally, it is clear that gravitational collapse in AdS (or confined grav- 
ity) is driven by some focusing mechanism by which the bulk energy density 
cascades to shorter scales. We have analyzed this behavior from both bulk 
and boundary perspectives. In the latter, we have uncovered an exponential 
behavior in the growth of the spectral bandwidth of the scalar field. This 
observation is also confirmed by examining the behavior at the origin and, 
interestingly, 11^ displays a behavior consistent with a Kolmogorov-like mech- 
anism. However, further investigations are required to put this observation 
on firmer footing. 

Acknowledgments : It is a pleasure to thank Oscar Dias, Chad Hanna, 
Gary Horowitz, Pavel Kovtun, Robert Myers and Jorge Santos for interest- 
ing and helpful discussions. This work was supported by the NSF (PHY- 
0969827 to Long Island University) and NSERC through Discovery Grants 
(to AB and LL) as well as CIFAR (to LL). LL and SLL thank the KITP for 
hospitality where parts of this work were completed. Research at Perimeter 
Institute is supported through Industry Canada and by the Province of On- 
tario through the Ministry of Research & Innovation. Computations were 
performed thanks to allocations at the Extreme Science and Engineering 
Discovery Environment (XSEDE), which is supported by National Science 
Foundation grant number OCI-1053575 as well as SHARCNET. 



References 

[1] P. Bizon and A. Rostworowski, "On weakly turbulent instability of 



anti-de Sitter space," 


Phys. Rev. Lett. 107 (2011) 031102| 


[arXiv: 1104.3702 [gr-qc] 





[2] J. Jalmuzna, A. Rostworowski, and P. Bizon, "A Comment on AdS 
collapse of a scalar field in higher dimensions," 
{Phys. Rev. D84 (2011) 085021, ,arXiv: 1108.4539 [gr-qcT 



[3] F. Pretorius and M. W. Choptuik, "Gravitational collapse in 



(2+l)-dimensional AdS space-time," \Phys.Rev. D62 (2000) 124012 
[arXiv:gr-qc/0007008 [gr-qc] 



[4] D. Garfinkle and L. A. Pando Zayas, "Rapid Thermalization in Field 
Theory from Gravitational Collapse," \Phys.Rev. D84 (2011) 066006] 
farXiv : 1 106 . 2339 [hep-th] . 



31 



[5] D. Garfinkle, L. A. Pando Zayas, and D. Reichmann, "On Field 
Theory Thermalization from Gravitational Collapse," 



JHEP 1202 (2012) 119, arXiv: 1110.5823 [hep-thj] 



[6] O. J. Dias, G. T. Horowitz, and J. E. Santos, "Gravitational Turbulent 
Instability of Anti-de Sitter Space," 

[Class. Quant. Gray. 29 (2012) 194002 , ^arXiv : 1 109 . 1825 [hep-th] 



[7] M. W. Choptuik, "Universality and scaling in gravitational collapse of 
a massless scalar field," fPhys. Rev. Lett. 70 (1993) 9-12^ 

[8] C. Gundlacli and J. M. Martin-Garcia, "Critical phenomena in 
gravitational collapse," Living Rev.Rel. 10 (2007) 5, 
,arXiv: 071 1.4620 [gr-qc] , 

[9] M. T. Anderson, "On the uniqueness and global dynamics of AdS 
spacetimes," \Class. Qu ant Gray. 23 (2006) 6935-6954^ 



farXiv:hep-th/0605293 [hep-th] 



[10] M. Dafermos. Private communication, (2011). 

[11] A. Ishibashi and K. Maeda, "Singularities in asymptotically anti-de 



Sitter spacetimes," 


Phys.Rev. D86 (2012) 104012 


arXiv: 1208. 1563 [hep-th] 





[12] O. J. Dias, G. T. Horowitz, D. Marolf, and J. E. Santos, "On the 
Nonlinear Stability of Asymptotically Anti-de Sitter Solutions," 
{Class. Quant. Gray. 29 (2012) 2350191 [arXiv : 1208 . 5772 [gr^gcT 



[13] A. Buchel, L. Lehner, and S. L. Liebling, "Boson stars in AdS," in 
preparation (2012) . 

[14] D. Astefanesei and E. Radu, "Boson stars with negative cosmological 



constant," Nucl.Phys. B665 (2003) 594-622 



arXiv:gr-qc/0309131 [gr-qc] 



[15] M. W. Choptuik, "Experiences with an adaptive mesh refinement 

algorithm in numerical relativity," in Frontiers in numerical relativity, 
C. Evans, L. Finn, and D. Hobill, eds. Cambridge University Press, 
1989. 

[16] L. Lehner, S. L. Liebling, and O. Reula, "AMR, stability and higher 
accuracy," ' Class. Quant.G ray. 23 (2006) S421-S446t 
farXiv:gr-qc/0510111 [gr-qc] [ 



32 



[17] M. Maliborski, "Instability of Minkowski spacetime enclosed in a 



cavity," 


Phys.Rev.Lett. 109 (2012) 221101, 


arXiv : 1208 . 2934 [gr-qc] 




http : //link. aps . org/doi/10 . 1103/PhysRevLett . 109 . 221101 



[18] M. Frigo and S. G. Johnson, "Fastest fourier transform in the west 
(fftw)." |http : //www ■ f f tw ■ org/l 2012. 



[19] H. de Oliveira, L. A. P. Zayas, and E. Rodrigues, "A 

Kolmogorov-Zakharov Spectrum in AdS Gravitational Collapse," 
[arXiv: 1209 .2369 [hep-th] | 



[20] A. Buchel, L. Lehner, and R. C. Myers, "Thermal quenches in N=2* 
plasmas," ^JHEP 1208 (2012) 049, ,arXiv: 1206 . 6785 [hep-th] ( 



33 



